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We present a framework for analyzing the rheology of dense driven granular materials, 
based on a recent proposal of a stress-based ensemble. In this ensemble fluctuations in a 
granular system near jamming are controlled by a temperature-like parameter, the angoric- 
ity, which is conjugate to the stress of the system. In this paper, we develop a model for 
slowly driven granular materials based on the stress ensemble and the idea of a landscape 
in stress space. The idea of an activated process driven by the angoricity has been shown 
by Behringer et al (2008) to describe the logarithmic strengthening of granular materials . 
Just as in the Soft Glassy Rheology (SGR) picture, our model represents the evolution of 
a small patch of granular material (a mesoscopic region) in a stress-based trap landscape. 
The angoricity plays the role of the fluctuation temperature in SGR. We determine (a) the 
constitutive equation, (b) the yield stress, and (c) the distribution of stress dissipated during 
granular shearing experiments, and compare these predictions to experiments of Hartley & 
Behringer (2003). 
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1. Introduction 

A striking feature of dry granular materials and other athermal systems is that they form 
force chain networks, in which large forces are distributed inhomogeneously into linear 
chain-like structures, in response to applied stress (Majmudar & Behringer (2005)). A 
number of recent experimental studies have visualized and quantified force chain networks 
using carbon paper (Jaeger et al (1996)) and photo-elastic techniques (Howell & Behringer 
(1999), Veje et al (1999)). These studies have demonstrated that the characteristics of force 
chain networks are acutely sensitive to the nature of the prepared state, especially near the 
jamming transition (Majmudar & Behringer (2005)). For example, in isotropically com- 
pressed systems, force chain networks are ramified with only short-ranged spatial cor- 
relations of the stress. In contrast, in sheared systems, aligned force chains give rise to 
long-ranged spatial correlations of the stress in the direction of the shear. Any theory of 
granular rheology has to, therefore, incorporate the effects of these grain-level structures 
on the macroscopic, collective response. 

It has been long realized that classical elasticity theory or Newtonian fluid dynamics 
is inadequate for describing the response of granular materials for a couple of important 
reasons (Bouchaud (2003), Cates et al (1998), Blumenfeld (2004)): (1) tensile stresses are 
completely absent in dry granular materials, and therefore the cohesion of the granular 
assembly is induced by the applied stress making the zero-stress state ill defined, (2) there 
is an indeterminacy of the forces at the microscopic level due to friction and disorder, (3) 
granular materials are athermal and dissipative, therefore there is no established statistical 
framework that bridges the structure at grain scales to a continuum elasto-plastic theory at 
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large length scales. In addition, granular solids often occur in isostatic states in which the 
number of degrees of freedom matches the number of constraints (Tkachanko & Witten 
(1999), Moukarzel (1998)). It has been shown that the critical properties of the isostatic 
point, not elasticity theory, determine the mechanical response of these marginal solids 
(Wyart (2005)). 

A framework, known as soft-glassy rheology has been used widely to understand the 
rheology of materials such as foams, colloids, and grains (Sollich (1998)). In this work, 
we use a recently proposed statistical framework, the stress ensemble for describing the 
response of static granular media (Henkes & Chakraborty (2009), Henkes et al (2007)), 
to formulate a theory of the response of slowly (quasistatic) sheared granular media and 
discuss its predictions for granular rheology. 

2. Stress Ensemble 

As shown previously (Henkes & Chakraborty (2009), Ball & Blumenfeld (2002)), the 
stress field of a a mechanically stable granular material can be fully described by a spatially- 
varying scalar field \|/ in two dimensions, and a tensorial field in higher dimensions (Henkes 
& Chakraborty (2009)). These fields can be used to establish a rigorous conservation prin- 
ciple, valid for grain assemblies in mechanical equilibrium, at the grain level in two di- 
mensions and at a continuum level in higher dimensions (Henkes & Chakraborty (2009)). 
This conservation principle, in conjunction with a maximum entropy hypothesis (Jaynes 
(1957)), leads to a generalization of equilibrium thermodynamics to the ensemble of me- 
chanically stable granular states (Henkes & Chakraborty (2009), Henkes et al (2007)) 
where the new conserved tensorial quantity, E to be discussed below, plays the role of 
energy. The microcanonical version of this ensemble is characterized by the complexity, 
a measure of the number of grain configurations compatible with a given value of E. The 
assumption of entropy maximization leads to the definition of a temperature-like intensive 
variable, the angoricity (Henkes & Chakraborty (2009), Blumenfeld & Edwards (2003)), 
which is a tensor. The predictions of this stress ensemble have been remarkably successful 
in describing results of simulations (Henkes et al (2007), Lois et al (2009)) and experiments 
in both static (Majmudar & Behringer (2007), Lois et al (2009)) and slowly driven gran- 
ular media (Behringer et al (2008)). As an example, the stress ensemble framework and 
the constraint that all contact forces are non-negative has been used to develop a Ginzburg- 
Landau type action for jammed granular systems in terms of \|/. The theory has been used 
to predict the spatial correlations of stress in systems subjected to isotropic compression 
(Henkes & Chakraborty (2009)), and shear (Lois et al (2009)), and to construct a mean- 
field theory of unjamming under isotropic compression (Henkes & Chakraborty (2009)). 
In slowly sheared granular systems, the observation of a logarithmic strengthening has 
been explained by the stress ensemble approach(Behringer et al (2008)). In this paper, we 
present a detailed development of the theory of slowly sheared granular matter, and make 
falsifiable predictions. Preliminary tests of the theory, based on comparison to experiments 
using photo-elastic disks (Hartley & Behringer (2003)) are also presented. 

(a) Conservation principle and Complexity 

To understand the mechanical response of granular materials, one needs a theoretical 
approach that can bridge the gap between microscopic, grain-level quantities and macro- 
scopic, collective properties. Fluctuations are inherently related to the number of micro- 
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scopic states available under a given set of macroscopic parameters. In equilibrium ther- 
modynamics, the microcanonical entropy, or its derivatives in other ensembles, is the mea- 
sure we use to calculate fluctuations and response. Conservation of energy allows for a 
rigorous definition of this measure since the states with different energies are not mixed 
by dynamics. In disordered systems such as spin glasses, the concept of complexity has 
been useful in formulating a framework for understanding collective properties (Bouchaud 
et al (1996)). Complexity is a measure of the number of states associated with a free en- 
ergy minimum, and has proven to be a useful concept for disordered systems that have 
many metastable states. In mean-field models, these minima are separated by barriers that 
diverge in the thermodynamic limit, and one can in principle count the number of states 
unambiguously. Can we identify a physical variable in granular materials, which also have 
many metastable states (Bouchaud (2003)), which is conserved by natural dynamics and, 
therefore, leads to a formalism for defining complexity? The lack of a natural dynamics in 
granular systems makes this a difficult proposition. In mechanically stable states, however, 
there is a topological conservation law that allows us to proceed to define the analog of 
complexity. The topological nature implies that a change in this physical variable can be 
achieved only through rearrangements that involve the boundaries or the whole system. 
A gedanken experiment serves to illustrate the conservation law. If we draw an imagi- 
nary line through a grain assembly and calculate the value of the total normal force being 
transmitted along this line, then force balance ensures that this value remains unchanged 
as this imaginary line is translated perpendicular to itself across the assembly (Metzger 
(2008)). A formal theoretical framework can be formulated using the force-moment ten- 
sor E = / d d rd(r), where 6(r) is the local stress tensor. It has been shown (Henkes & 
Chakraborty (2009), Henkes et al (2007)) that E depends only on the boundary conditions 
of the packing. Thus, the phase space of all mechanically stable configurations can be di- 
vided into sectors labeled by E. Configurations in different sectors are disconnected under 
any local dynamics, and E plays a role similar to total energy in equilibrium thermodynam- 
ics or free-energy minima in spin glasses. Based on some very general principles, such as 
the factorization of states (Bertin et al (2006)), the conservation of E has been shown to 

lead to an intensive variable, a = -|^p, where 5(E) is the entropy of the sector labeled by E 
(Henkes & Chakraborty (2009)). The inverse of a, the angoricity (Edwards & Blumenfeld 
(2007)) is expected to play the same role as temperature in equilibrium thermodynamics, 
if processes that create granular assemblies achieve entropy maximization. Comparisons 
with simulations have demonstrated that a mechanically stable assembly of grains has the 
same value of angoricity throughout its interior, while 6(r), the local stress fluctuates. 
The probability of occurrence of a microscopic state v within a mechanically stable grain 
packing, maintained at an angoricity, a, is: 

P v = (l/Z)e- &± v , (2.1) 

analogous to the Boltzmann distribution (Henkes & Chakraborty (2009)). The state v = 
{ri,fjj}, where rj denote the positions of the grains, and fy the set of contact forces. 

A natural question to ask is whether the similarities between temperature and angoricity 
extend to dynamics of slowly deformed granular systems, situations analogous to dynam- 
ics of thermal systems close to thermal equilibrium where stochastic dynamics such as 
relaxational, Langevin dynamics provides a good description (Goldenfeld (1992)). In the 
granular context, we are restricting our attention to systems that are close to granular equi- 
librium defined by Eq. 2. 1 . As in systems close to thermal equilibrium, we can construct a 
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stochastic, Markov process based on the concept of detailed balance: 

= * (2.2) 

This is an additional assumption beyond the ones entering the construction of the stress 
ensemble that leads to the distribution of mechanically stable states, P v . Stochastic equa- 
tions, based on detailed-balance conditions lead to the Boltzmann-like distribution, Eq. 
2.1, as the time-independent distribution of states (Goldenfeld (1992)). If, however, the 
stress-landscape (or energy landscape in a thermal ensemble) is such that complete equili- 
bration is not possible then the stochastic equations lead to a rate of escape over a barrier 
AE, which is the well-known Kramers rate: g~" :Al: (Kramers (1940)). A simple example 
of a landscape where equilibration is not possible is one that has a single local minimum 
and a single local maximum. In that situation the system "escapes" over the barrier with 
a rate prescribed by Kramers but cannot equilibrate. As will become clear from the dis- 
cussion below, trap models describe systems that escape over barriers with the Kramers 
rate, however, subsequent sampling of traps is described by a quenched distribution and 
not by detailed balance. We adopt this form of activated dynamics to analyze the rheology 
of dense, quasistatic granular flows. 

Our framework of granular rheology is, therefore, based on a rigorous conservation 
principle and postulates, similar to the ones adopted in the development of equilibrium 
thermodynamics (Callen (1957)). These postulates are (a) factoris ability of the dynamics- 
dependent frequency with which each grain configuration is accessed (Henkes & Chakraborty 
(2009)), which allows for the definition of the angoricity, (b) maximization of entropy, 
which implies equality of the angoricity across a system in granular equilibrium and a 
Boltzmann-like distribution (Eq. 2.1) (Henkes & Chakraborty (2009)), and (c) satisfaction 
of detailed balance by the rates of transition between microscopic states. 



3. Trap models and rheology 

As discussed earlier, the existence of a large number of different microscopic metastable 
states that are macroscopically equivalent puts granular materials into a wider class of 
previously well studied systems, including gels, glasses, colloids, emulsions, polymers 
and foams. A minimal model that encapsulates these features, and leads to a glass tran- 
sition is the "trap model" (Monthus & Bouchaud (1996)). In this model, one considers a 
system made of independent subsystems of a certain size where each subsystem acts 
self-coherently and independent of the others. Their dynamics involve hopping between 
different metastable states aided by some kind of fluctuation. This idea of a random walk 
in a rugged landscape has its roots in the context of glasses (Bouchaud & Georges (1990)). 
Below, we discuss the trap model and its generalization, which is the framework of soft 
glassy rheology (SGR) (Sollich (1998)) in some detail. 



(a) Bouchaud 's trap model for glasses 

Monthus & Bouchaud (1996) constructed a one-element model for glasses. In this 
model, there exists an energy landscape of traps with various depths E. An element hops 
between traps when activated, where the fluctuation is assumed to be thermal. At a temper- 
ature kgT = 1 /P, the probability of being in a trap with depth E at time t evolves according 



Article submitted to Royal Society 



Dynamics in a Stress Landscape 



5 



to: 

^-P{E,t) = -me-V E P{E,t) + (S){t)p{E) (3.1) 

where, 

co(f) = COo (e~ p£ ) P (E,t) = wo j dEP(E,t)e-V E (3.2) 

The first term on the right-hand-side (rhs) of Eq. 3.1 is the rate of hopping out of a trap, 
where e~$ E is the activation factor and COo a frequency constant. The model assumes that 
choosing a new trap is independent of history, so that a new trap is chosen from a distribu- 
tion of trap depths p(E) that reflects the underlying disorder in the glassy models. Eq. 3.2 
gives the average hopping rate. This rate, multiplied by the distribution of trap depth, Eq. 
3.1 gives the probability rate of choosing a new trap. 

The existence of a glass transition in this model can be demonstrated as follows: assume 
the distribution function has a simple exponential tail such as: 

p(E) - e-^ E , (3.3) 

where Po is a fixed parameter describing the disorder in the inherent energy landscape. The 
physical justification for having an exponentially decaying p(E) is borrowed from systems 
with quenched random disorder, such as spin glasses, which use extreme value statistics 
(Monthus & Bouchaud (1996)) (We will later justify applying the same form of distribu- 
tion to granular materials). Then we can easily solve for the steady-state (^- t P(E,t) = 0) 
solution: 

P eq (E) - e p£ p(£) - e (Po-P)£ (3 .4) 

Immediately, we see that Eq. 3.4 is non-normalizable for p > po or T < Tq. This shows 
that below a temperature 7b, the system is out of equilibrium; it is non-ergodic and ages by 
evolving into deeper and deeper traps. Therefore, we call T — To a point of glass transition. 
(Monthus & Bouchaud (1996)). 

It should be noted that dynamical equations of the trap model do not lead to a Boltz- 
mann distribution at T since the equations do not obey detailed balance, and specifically, 
the traps are sampled from a quenched distribution. The escape rate from a trap is, however, 
determined by the Kramers process. 

(b) Soft Glassy Rheology 

In (Sollich (1998)), a simple phenomenological model termed "Soft Glassy Rheology" 
(or SGR) was proposed to describe the anomalous rheological properties of "soft glassy 
materials". The model incorporated Bouchaud's glass trap model and an extra degree of 
freedom, the strain. The model assumes a macroscopic soft glassy material can be sub- 
divided into a large number of mesoscopic regions, each having a linear size With each 
mesoscopic region one can then associate an energy E and a strain variable /, which both 
evolve with time. The number of mesoscopic regions must be made large enough so that 
ensemble averaging can be performed over them to yield macroscopic properties of soft 
glassy materials. Similar to Bouchaud's trap model, the SGR model assumes the meso- 
scopic regions "live" in an energy landscape. This is a mean-field energy landscape in that 
it is not characterized by a metric, and is defined only by the distribution p(E). The new 
strain variable describes local elastic deformation of the mesoscopic regions, so that / con- 
tributes quadratically to the energy E. Since a strain variable was added to the SGR, the 
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model can now describe material under imposed shear strain, also in a mean-field spirit, all 
mesoscopic regions respond uniformly to externally imposed shear strain. Different from 
Bouchaud's trap model, thermal fluctuations are unimportant in soft materials (IcbT is too 
small to cause structural rearrangements). In the SGR model, it is the fluctuations in the 
elastic energy that facilitate structural rearrangements. This fluctuation is determined by 
a temperature-like quantity x called the "noise level". As a result, the escape rate from a 
trap becomes e~ E l x . The use of x is a mean-field approach to describing all interactions 
between mesoscopic regions. 



(c) Incorporating the stress ensemble into the SGR framework 

The idea of using a temperature-like quantity to replace thermal fluctuations in the 
SGR model is reminiscent of the stress ensemble for granular materials, where the inverse 
angoricity, 1 /a, plays the role of the noise, x, in SGR. The activated process in the stress 
ensemble is a result of a coupling between different mesoscopic regions through stress 
fluctuations. Each mesoscopic region can be viewed as existing in a bath of stress fluctua- 
tions, which is characterized by the angoricity, and these fluctuations can lead to activated 
processes analogous to those occurring in thermal systems. Since the granular material in 
relevant experiments (Hartley& Behringer (2003)) are sheared setups, and to simplify the 
model, we will use a scalar model that incorporates only the shear components of the stress 
through r, the deviatoric part of £, and a, the shear component of the inverse-angoricity 
tensor (referred to as shear-angoricity below), which represents the bath of stress fluctua- 
tions that the mesoscopic region is in contact with. 

To adopt the SGR framework to the stress ensemble, we need to define metastability 
and escape processes in stress rather than energy space. Since the stress ensemble is based 
on the premise that there are many states {r^fy} with a given F, and that these states are 
broadly sampled even if not with equal weights (Henkes, 2009), we assume that mesco- 
scopic regions of granular assemblies can be in metastable equilibrium characterized by 
the Boltzmann-like distribution (Eq. 2.1). A dynamics obeying detailed balance (Eq. 2.2) 
would lead to this distribution. As mentioned earlier, the trap model and SGR, however, 
describe systems in the presence of quenched disorder which define the distribution barrier 
heights (or equivalently trap depths) Ar that can be crossed in a finite time. The system 
can escape from a metastable state represented by a trap. Stochastic equations based on the 
stress ensemble allows us to adopt the Kramers' escape rate approach. Replacing temper- 
ature by the shear angoricity, the escape rate is therefore, e~ aAr , which replaces e~$ E of 
Bouchaud's trap model. As in the original trap model framework, the exploration of the 
traps in stress space, following this escape is assumed to be determined by the intrinsic, 
quenched-in distribution of trap depths, and this dynamics does not lead to the Boltzmann- 
like distribution, Eq. 2.1 in the long-time limit. The quenched-in distribution of trap depths 
determines whether complete equilibrium is attainable in the absence of shearing (Monthus 
& Bouchaud (1996)). The activated process in the original trap model was thermal, and in 
SGR the activation is controlled by the noise x. In the scalar version of the stress-ensemble, 
the activation is in stress space and is controlled by the angoricity a. The disorder in grain 
packings is represented by an intrinsic distribution of trap depths, just as in SGR or the 
original trap model. 

We can now write a generalization of the SGR model to granular materials using the 
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stress ensemble: 

^p(r m ,r,o = -co oe - a(r '"- r) p(r m ,r,f) + co(f)p(r ffl )8(r) (3.5) 

at 

Eq. 3.5 describes the evolution of the probability of finding a mesoscopic region with a 
"yield stress" F m (again, strictly F m has the unit of force moment) and an instantaneous 
stress r. The yield stress F m is used to define the depth of the trap the system is in, and 
since F is the instantaneous stress of the system, the barrier is Ar = F m — F. The trap depth 
E of SGR or the trap model is therefore replaced in the stress ensemble by F m . 

The model accounts for the successive stress buildup and yielding events inside a meso- 
scopic region. When the instantaneous stress (due to external driving) is less than the yield 
stress, or F < F m , the mesoscopic region is considered to be inside a "trap" of depth F m . It 
is possible to yield below the yield stress because of the activated nature of the process. At 
this point the mesoscopic region "sees" a reduced trap depth (i.e. effective barrier height): 
F m — F, and hence can undergo an activated process to "hop" out of the trap with a rate pro- 
portional to e~ a ( r "~ r ). Another process through which the mesoscopic regions can yield 
is that if no activated process takes place while in the trap (r < F m ), F can increase due to 
external applied shear up to the point F = F m where the mesoscopic region yields, hence 
leaving the trap. This is the only possible process at the limit of zero angoricity (1/oc — ► 0), 
and is not explicitly included in Eq. 3.5. Either way, after the mesoscopic region leaves a 
trap, it choose a new trap according the distribution p(F m ). Fig. 1 shows a schematic of 
this dynamics of the mesoscopic region in the landscape. 

The mesoscopic regions are assumed to be large enough to be treated as deforming 
elastically until a "plastic" yield event occurs. In the trap picture, these plastic events are 
equated with the mescoscopic region leaving a trap. The grains are driven by a constant 
shear rate in the experiments (Hartley & Behringer (2003)). In a mean-field approach simi- 
lar to SGR, we assume this applied shear is transmitted without decay throughout the entire 
material. Thus each mesoscopic region is driven by the same shear rate. Since the meso- 
scopic regions deform elastically, the stress F increases linearly with time up to the point 
of yielding: 

o=lr = kl = kft (3.6) 

where o is the shear stress, S is the area of the mesoscopic region, and k is the elastic 
constant relating the shear stress to the shear strain, /. The continual process of stress 
buildup and yielding creates a "sawtooth" pattern in time as show in the bottom of Fig. 1. 
With a linear dependence of F on time, between yield events, we can rewrite Eq. 3.5: 

^p{F m ,F,t) = -5^p(r m ,r,f)-co oe - a(r '»- r) p(r m ,r,f) + co(f)p(r m )5(r) (3.7) 

The full time-differential on the l.h.s of Eq. 3.5 has been converted to a partial time deriva- 
tive plus an advective term. The delta function in the last term on the r.h.s. of Eq. 3.7 means 
that immediatly after choosing a new trap, the mesoscopic region is in a stress-free state. 
See points "3" and "5" in Fig. 1. If y = 0, the model reduces to Bouchaud's trap model in 
stress space. The steady state distribution, P(r m ,r) °< e~ a ( r '"~ r )p(r m )8(r), which is the 
equilibrium probability of finding a region with yield stress F m . Since there is no imposed 
shearing, the instantaneous stress is zero. 

Physically, hopping out of a trap equates to the grains in the mesoscopic region mak- 
ing a re-arrangement. These events are visible in experiments and if analysed with photo- 
elastic techniques, are just buckling or the breaking of force chains. After this re-arrangement, 
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r(t) 



t 

Figure 1. Schematic showing dynamics of the model in the stress landscape. The black dots with the 
numbering (1 thru 6) represent the state of the mesoscopic region at different times. In the top figure, 
each "trap" in the landscape is given by its depth T m . The amount of stress buildup, while still in the 
trap, is indicated by the height T from the bottom of the trap. For "1", the mesoscopic region has just 
fallen into the trap and it has no stress build up, so it has value of T = in the bottom figure. As 
constant shear stress is applied to the mesoscopic region, its stress starts to build up. In the top figure, 
"2" is at a height from the bottom of the trap; and in the bottom figure, "2" has increased linearly in 
value from "1". Next, via the activated process, the mesoscopic region makes a jump from one trap 
to another. Note the traps shown are intentionally made disjoint, this is to indicate that every trap can 
be accessed from every other trap and there is no specific connectivity. 

new force chains form. The newly formed force chains are characterized by a new yield 
stress r m . 

The distribution of trap depths p(F m ) describes the disorder in the stress landscape. 
Theoretically, this distribution can be measured in experiments where no activated stress 
yielding occurs, or when stress fluctuations exist only locally which leads to the angoric- 
ity (1/oc) being very small. This is the case when no neighboring stress fluctuation is 
"felt" within the mesoscopic region. This is realized experimentally (Peidong & Behringer 
(2009)). In this experiment, shear is applied to a small local region which can be considered 
a single mesoscopic region. The stress yielding is due entirely to the driving. The distribu- 
tion of yield stresses in this case is exactly p(F m ). Generally, in sheared granular systems 
(e.g. Geng & Behringer (2005)), we can get an idea for the form of p(F m ) by looking at 
the tail of the stress distribution (large F). This is because large F values are accessed only 
when an activated process has not occurred to make the the system yield. They are also 
states that sample large F m values. 

Activated dynamics in stress space is not an entirely new concept. It was proposed by 
Eyring (Glasstone et al (1941) ), and has been applied to analyze the velocity profile in 
dense granular flows (Pouliquen & Gutfraind (1996)). In our work, the stress-activation is 
a natural consequence of the ensemble based on angoricity. 

We can write down a general solution to Eq. 3.7 regardless of the distribution of trap 
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depth p(r,„). Eq. 3.7 can be simplified with a change of variable: 8r = F(t) — Skyt, which 
does not explicitly depend on time (away from the points of yielding). This eliminates the 
advective term in Eq. 3.7 which then becomes: 

|p(r m ,sr,o = -coo e - a(r " 1 - (8r+ ^ )) J p(r m ,8r,f) + (o(f) P (r m )8(8r+5^) (3.8) 

A better intuition can be gained by solving for 3.8 while ignoring the second term on its 
r.h.s., the solution is: 

p(r B ,sr,o =Po(r m ,8r) exp (-J^W( 5r + s **')) (3 . 9 ) 

The solution is just an exponential decay with the time interval replaced by a time integral. 
This is the "effective time interval" defined in (Sollich (1998)): 

Z(f/;8r)= f dt"e< sr+sv < t ") (3.10) 
Jt' 

In essence the behavior of an element in a trap is just that of a decaying process, <~ e^l x . 
For us, the time is replaced by the effective time interval Z(f , 0; 8r) that grows much faster 
than linearly and the mean lifetime is x = (Oo~ l e ar,n . After the stress collapse of the first 
buildup occurring at t' , the element escapes into another trap chosen from the distribution 
p(r m ) and starts to undergo its own decay process via the effective time. Therefore the full 
solution is: 

p (r m , sr, t ) =p Q (r m , sr) exp (- ^z(t,o-, sr)) p.i i) 

+ £ dt'(o{t')p (r m ) 8 (sr + skjt') exp (- J^z (r , t '■ sr) ) 



4. Theory predictions and comparison to experiments 

Slow shearing experiments correspond to the long-time or steady state limit of our model. 
In such a limit we find a solution that is time independent: 

P(r m ,r) = p(r ffl )exp (-^f^ (e aT - 1)) , (4.1) 

subject to normalization. We now have a ensemble distribution function of Y with which 
we can calculate the constitutive relation by taking the ensemble average of T. Using Eq. 
3.6 we can write: 

o= l -{T) P[rmX) = l -j~drrf~dr m p{r m ,r) (4.2) 

We first consider the simplest case of a distribution of trap depths, p(r„,) = 8(T m — To) 
with To > 0. Although not a physical description, this choice will prove to be a fruitful 
exercise. Eq. 4.2 becomes: 

1 jVdrr -Umru- " 
aS E X {A) { ' 
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where the constant A = jpf^e ar ° and E\ is the exponential integral with n = 1. In the 

interesting limit of A — > or y » coo e~ ar °/ Ska, we can expand Eq. 4.3 to 0(A°) and 
obtain: 

° = 2o5 l0g (Y) + 2o5 ^ A°° ) + aF ° ~ 2 Yc) ( 4 - 4 ) 
where y c is the Euler-Mascheroni constant ~ 0.577. The leading behavior of the consti- 
tutive equation has a logarithmic dependence on the shear rate. This agrees qualitatively 
with experiments (Behringer et al (2008), Hartley & Behringer (2003)). In the opposite 
limit A — > oo or y « coo e~ ar °/ Ska, we can expand Eq. 4.3 to 0(1 /A) and obtain: 

ke ar 

a = y (4.5) 

This reveals a Newtonian regime with viscosity r| = k e aF °/(Oo. 

(a) Constitutive equation for the exponentially decaying distribution of trap depths 

As discussed in the previous section, an exponentially decaying p(T m ) is a better de- 
scription of granular experiments as well as spin-glasses. We use: 

p(r B ) = -U- a ° r '" (4.6) 

a 

where ao is a constant measuring the disorder in the stress landscape. Inserting Eq. 4.6 into 
Eq. 4.2 we obtain: 

, r°° dW \og(W+(x>o/Skya) , w s 
a = — log(Skya (£>o) H — mr 1 — — (4-7) 

°" a ^ Jo ww+^o/sm y{ - x ' w > 

where x is given by the ratio: x — Oo/a, and y(x, W) is the lower incomplete gamma func- 
tion. Again we look in the limit of y » a>o Ika and the constitutive equation is approxi- 
mated to zeroth order as: 

This shows that the logarithmic dependence of shear stress on shear rate is independent 
of the form of the distribution p(r m ) used. In Fig. 2 we compare this result with the data 
obtained in 2D Couette granular shearing experiments (Hartley & Behringer (2006)). 
In the opposite limit, y « a>o I Ska we again obtain a Newtonian regime: 

k x—1 . , , „ N 

o = -y (4.9) 

(Oqx — 2 

The viscosity diverges atx = 2, the same result is obtained in SGR (Sollich (1998)). 

(b) Yield Stress 

When y— 0, the stress has a finite value for ao < a or x < 1. This is the yield stress, 
and it vanishes at the same point as the glassy transition point in Bouchaud's trap model. 
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0.01 0.1 1 10 100 

Shear rate, y 

Figure 2. Linear-Log Plot of mean stress in a segment of 2D granular Couette experiment contain- 
ing roughly 200 particles vs. shear rate. Data, taken from (Hartley & Behringer (2003)), are given 
for different packing fractions, which are given relative to the critical packing fraction (|) c where 
the system loses all rigidity. Various symbols are: upward triangles: = <|) c + 0.00035; diamonds: 
<|> = <j> c + 0.00141; squares: § = (|) c + 0.00528; solid disks: = ty c + 0.00528; downward triangles: 
<|> = $ c + 0.00800; and hollow circles: <|) = <|> c + 0.01373. 



The yield stress can be calculated by inserting y = into Eq. 4.2. Performing the T integral 
first in Eq. 4.2, we obtain: 

. r^/Skja, x -l G() 

o= V° - dZZ ° [Z) , (4.10) 
aSf*o/s k ya dzzX _ lEi{z y 

where G(z) is the Meijer G-function: 

°u o!o! )=r^ ex p(-^- 1 )) 

When y= 0, Eq. 4.10 can be calculated exactly: 



a y = a(y = 0) = -kxy c -kx— log (r(x) ) , (4.11) 



A plot of the yield stress is shown in Fig. 3 



(c) The distribution of stress drops 

It is often easier to measure distribution of stress drops in experiments than the stresses 
themselves, and the distribution provides a more stringent test of theoretical frameworks. 
While in a trap with depth F m , the probability to build up the stress to a value of F is given 
by: (c/Eq. 4.1 ) 

P S (T) - exp (- - &0 nr (e ar -l)) 
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Figure 3. Theoretical Plot of the yield stress Eq. 4.1 1 as a function of x = do/a. It vanishes at x=l. 

On the other hand, as the stress grows, it is increasingly likely to fail due to the activated 
process. The rate of failure is given by: 

Pf(r) - e - a ( r »'-n 

Then, the probability of building the stress up to F and then to fail at this point is given 
by the product P s x Pf. Since the stress falls to zero after a failure event, we can call the 
magnitude of the stress drop AT\ Its distribution for a given p(T m ) is proportional to: 

P(AT) oc P s (Ar)P f {AT) 

Also sum over all possible traps: 

P(AT)« Jdr m p{r m )P s {AT)P f {Ar) (4.12) 

Using the exponentially decaying distribution p(T m ) = e~ a o r ™, we obtain the normalized 
distribution: 

where we see the lower incomplete gamma function again. There are three regimes in 
which Eq. 4.13 can be simplified. First, using the limit y{$,y) — ► 1/P asy — > 0, we get: 

P(AT) = (^-) £L e aAr when ^( e aAr -l)^0 (4.14) 
v ; \x+lJSkj Skja y ' 

The opposite limit for the lower incomplete gamma function is y(P,>') — > (P — 1) ! as y — > °o 
for which we get: 

(4.15) 
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Table 1. Results of fitting Eq. 4.15 to the experimental data of Hartley & Behringer (2003) 
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Finally, the more extreme case of Eq. 4.15 as Ar — > °°: 

P(Ar) = au! f-^-) e- a ° Ar w/!en Ar -» °° (4.16) 
v /c 'yet y 

In Eq. 4.16, we recover the exponential distribution of depths. 

With an analytical form of P(Ar), we can compare and fit to data obtained by Hartley & 
Behringer (2003). The data is fitted to the form Eq. 4.15 and therefore contains two fitting 
parameters: a and oco- Table 1 shows the fitting results. In Fig. 4, four fitting results are 
plotted and compared with data. As seen from Table 1, both a and oco are functions of the 
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(e) (d) 

Figure 4. Log-Log plot of experimental distribution of stress drops and its fitting to Eq. 4.15. The 
experimental and fitting parameters are: (a) y = 0.3318ra/fz, <|> — (|) c = 0.0053, a = 0.70, do = 0.3247, 
(b) y = 0.9972m//z, §-§ c = 0.0053, a = 0.7, cc = 0.2660 , (c) y= 53155mHz, (|> - ty c = 0.0088, 
a = 1.2, a = 0.2131, and (d) y = 0.0657mHz, ty-ty c = 0.0137, a = 0.4, a = 0.2547. 

shear rate and the packing fraction. The data fitting gives us insights into the relationship 
between the shear-angoricity (x = GCn/ a), the shear rate (y), and the packing fraction <|> — <j) c . 
We can see this relation in a plot of x vs. y in Fig. 5 for different (j) — (j) c . The first observation 
one makes is that x is controlled by both y, and — (j) c , implying that the SGR assumption 
of x being independent of the shear rate does not apply to these measurements in sheared, 
granular media. In addition, over a decade of (j) -(j) c values, and three decades of y, x seems 
to be described by a scaling form: 




(4.17) 



with A ~ —0.4, and the scaling function (for (j) > (j) e ) f+(z) being a decreasing function 
of its argument, z. The dependence on y is, however, weaker than the dependence on (j) — 
(j) c . Although better statistics are needed to pin down the scaling form, it is intriguing to 
examine its consequences as such. The scaling itself hints at ((j) — (j) c = 0,y = 0) being a 
critical point, which is consistent with the observations of Olsson & Teitel (2007), and 
the idea of Point J being a special point in the jamming phase diagram (Liu & Nagel 
(1998)). At zero shear rate, x becomes independent of (j) — (j) t , and the scaling form is, 
therefore, consistent with the packings possessing a yield stress for all (j) > (j) t . From the 
perspective of x being the fluctuation temperature, the observation is that stress fluctuations 
grow as (]) approaches (j) c from above. Experiments do indicate growing fluctuations as a 
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system approaches the unjamming transition below which it cannot sustain shear (Howell 
&Behringer(1999). 
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Figure 5. Log-log plot of the scaled the relationship (Eq. 4.17). The scaling exponent is found to 
be: A = —0.4 The data points are: disks: — (|) c = 0.0004, squares: — <|) c = 0.0014, diamonds: 
<t> — 4*c — 0.0053, upward triangles: — (|) c = 0.0053, and downward triangles: — <|) c = 0.0053. 
Inset: Unsealed log-log plot of x vs. j. 



5. Conclusion 

A stress-based statistical ensemble has been used in conjunction with the idea of a stress 
landscape to construct a model for granular rheology. Comparisons of the model to exper- 
iments on slowly sheared granular media indicate that the theoretical framework provides 
a semiquantitative description of the stress-response. The framework is an extension of 
the Soft Glassy Rheology approach, and the stress ensemble provides a natural definition 
for the fluctuation temperature. This fluctuation temperature can be obtained by fitting the 
theoretical predictions to experiments, and the results show that fluctuations are largest at 
the special point, Point J, where the shear rate goes to zero and the packing fraction ap- 
proaches the critical value below which one cannot construct a mechanically stable state 
of grains. In future work, we intend to extend the model to include volume fluctuations 
since Reynold's dilatancy is one of the signature properties of granular matter (Reynold 
(1885)). A fundamental premise of the stress ensemble framework is that just as in thermal 
equilibrium temperature is equalized, for a system in granular equilibrium, the angoricity 
is the same everywhere inside the system. Comparison to simulations has shown this to be 
true for a system under pure compression (Henkes et al (2007)). If we assume the same 
to be true for the shear-angoricity x, then according to Eq. 4. 17, the local packing fraction 
has to adjust to the local shear rate and regions of increased shear rate have lower packing 
fractions. An interesting question to ask is whether this is the origin of Reynold's dilatancy. 
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